Magnetoconductance of interacting electrons in quantum wires: Spin density 

functional theory study 
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We present systematic quantitative description of the magnetoconductance of split-gate quantum 
wires focusing on formation and evolution of the odd (spin-resolved) conductance plateaus. We 
start from the case of spinless electrons where the calculated magnetoconductance in the Hartree 
approximation shows the plateaus quantized in units of 2e^ /h separated by transition regions whose 
width grows as the magnetic field is increased. We show that the transition regions are related to the 
formation of the compressible strips in the middle of the wire occupied by electrons belonging to the 
highest (spin-degenerate) subband. Accounting for the exchange and correlation interactions within 
the spin density functional theory (DFT) leads to lifting of the spin degeneracy and formation of the 
spin-resolved plateaus at odd values of /h. The most striking feature of the magnetoconductance 
is that the width of the odd conductance steps in the spin DFT calculations is equal to the width 
of the transition intervals between the conductance steps in the Hartree calculations. A detailed 
analysis of the evolution of the Hartree and the spin-DFT subband structure provides an explanation 
of this finding. Our calculations also reveal the effect of the collapse of the odd conductance plateaus 
for lower fields. We attribute this effect to the reduced screening efficiency in the confined (wire) 
geometry when the width of the compressible strip in the center becomes much smaller than the 
extend of the wave function. A detailed comparison to the experimental data demonstrates that the 
spin-DFT calculations reproduce not only qualitatively, but rather quantitatively all the features 
observed in the experiment. This includes the dependence of the width of the odd and even plateaus 
on the magnetic field as well as the estimation of the subband index corresponding to the last resolved 
odd plateau in the magnetoconductance. 

PACS numbers: 73.21. Hb, 73.43.Qt, 73.43.Cd, 73. 23. Ad 



I. INTRODUCTION 

The quantized conductance of a two-dimensional elec- 
tron gas (2DEG) in the quantum Hall regime has gener- 
ated a tremendous attention since its discovery in 1980^. 
For a theoretical description of the integer quantum Hall 
(IQH) effect, the concept of edge states combined with 
the Landauer-Buttiker formalism is widely used^ii^. This 
approach is proved to be especially appealing for the 
description of electron transport in the quantum Hall 
regime in confined geometries like quantum wires or 
quantum point contacts (QPCs)^. 

Some aspects of the quantized conductance in the con- 
fined geometries in the IQH regime can be understood in 
a one-electron picture. This includes, for example, mag- 
netic depopulation of the subbands in a quantum wire^, 
and selective population and detection of edge channels 
by QPCs resulting in the observation of anomalous IQH 
effeclj^i^. In the one-electron picture the two-terminal 
magnetoconductance G of a quantum wire or a QPC ex- 
hibits quantized plateaus in units of 2e^ /h (for spinless 
electrons) separated by transition regions of an essen- 
tially zero width. The experiments, however, show that 
an extend of these transition regions can be comparable 
to the width of the plateaus^i^i^ii. This indicates that an 
accurate description of the magnetoconductance in the 
IQH regime even without accounting for spin effects re- 
quires approaches that go beyond a simple one-electron 
picture of non-interacting electrons. A quantitative elec- 



trostatic theory of interacting electrons in quantum wires 
was proposed by Chklovskii, Matveev and Shklovskii^. 
They demonstrated that in a strong magnetic field, al- 
ternating strips of compressible and incompressible liq- 
uids are formed in the center of the wire. They also 
evaluated the two-terminal magnetoconductance of the 
wire. In contrast to the one-electron description, the 
magnetoconductance of interacting electron was shown 
to exhibits very narrow quantized plateaus separated by 
much broader rises where the conductance was not quan- 
tized. This conclusion (being opposite to the prediction 
of the one-electron picture) is also in apparent disagree- 
ment with the experiments. This indicates that even for 
spinless electrons in the IQH regime an accurate quan- 
titative description of the magnetoconductance requires 
many-body quantum mechanical treatment. 

At low temperature in clean high-mobility samples 
the spin degeneracy is lifted and the additional plateaus 
at odd values of e'^/h become resolved^ii. This is 
due to the exchange and correlations effects leading 
to the strong enhancement of the electron g factor 
above its bulk valuo^. The effect of the many-body 
interactions on the spin splitting in quantum wires 
in the IQH regime has been a subject of numer- 
ous Studie9iaaiii2a3i4^i6^a8^^^^^. Thesc stud- 
ies have focused on various aspects of 2DEG in con- 
fined geometries including the structure of compress- 
ible/incompressible strips, suppression or enhancement 
of the g factor, subband spin splitting, spatial spin sepa- 
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ration and other. We however have not been able to find 
in the Hterature any systematic quantitative treatment of 
the magnetoconductance of the structures at hand. Sur- 
prisingly enough, even after two decades of the studies 
of IQH systems in the confined geometry, the question of 
the fundamental importance addressing the formation of 
the odd plateaus in the magnetoconductance and corre- 
sponding quantitative description of the plateau widths 
still remains unanswered. As discussed above, the struc- 
ture of the magnetoconductance plateaus remains poorly 
understood even for the case of spinless electrons. Recent 
advances in the field such as demonstration of the Mach- 
Zehnder— , Aharonov-Bohm^^ , and Laughlin quasiparti- 
cle interferometers^^ or prospects of the topological quan- 
tum computing^ has led to a renewed interested to the 
magnetoconductance in the quantum wires and related 
structures. Even though many of the the above systems 
operate in the fractional quantum Hall regime where the 
correlation effects become dominant, a detailed under- 
standing of the magnetoconductance in the IQH regime 
is the necessary prerequisite for the understanding of the 
magnetotransport in the fractional regime. 

In our previous publications we provided a systematic 
quantitative description of the structure and spin polar- 
ization of edge states and magnetosubband evolution in 
the quantum wire based on the self-consistent Green's 
function techniques combined with the spin density func- 
tional theory (DFT) ^^i^°'^^ or Hartree-Fock approacb^^. 
The main aim of this paper is to present a systematic 
quantitative description of the two-terminal magneton- 
ductance of the quantum wire with the focus on the for- 
mation and evolution of the exchange-induced odd con- 
ductance plateaus. The motivation for the present paper 
is the recent experimental studies of the spin-resolved 
magnetoconductance of the narrow channels in the IQH 
regime^. One of the remarkable finding of this experi- 
ment is the collapse of the spin splitting in the confined 
geometries for lower field. The spin-DFT magnetotrans- 
port calculations presented in this paper not only capture 
essential features observed in the experiment, but demon- 
strate rather good quantitative agreement with the calcu- 
lated and observed magnetoconductances. This includes 
the width and the position of the magnetoconductance 
plateaus (both odd and even), as well as predictions for 
the critical magnetic field where the odd plateaus dis- 
appear in the magnetoconductance. We therefore con- 
clude that the spin DFT approach represents the pow- 
erful tool to study large realistic quantum Hall systems 
containing hundreds or thousands of electrons, providing 
detailed and reliable microscopic information on wave- 
functions, electron densities and currents as well as the 
conductance. 



II. BASICS 

We consider an infinitely long wire in a perpendicular 
magnetic field B, see Fig. [TJ The bare electrostatic con- 
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FIG. 1: (a) Schematic diagram of a split-gate quantum 
wire, (b) The bare electrostatic parabolic confining potential 
Vconfiv) (Eq. Q) with Vo = -0.4 eV, hwo = 3.8 meV (solid 
lines). Dashed lines shows the calculated electrostatic con- 
finement for a realistic split-gate quantum wires depicted in 
(a) with the distance between the gates a — 500 nm, the gate 
voltage Vg — —2 V, 5— donor concentration = 6 ■ 10^'' m~^, 
the donor and electron distance to the surface are 100 nm and 
200 nm respectively and the Schottky barrier Vschottky ~ 0.8 
eV (The potentials from the gates and the donors are calcu- 
lated on the basis of Eqs. (2) and (3) of Ref. [l9h . 

finement (due to the split gates, the donor layer and the 
Schottky barrier) can be approximated by a parabolic 
potential, 

V,ons{y)^Va + — {u^y)\ (1) 

where Vq is the bottom of the potential, loq defines the po- 
tential slope, and m* = O.OGTme is the electron effective 
mass in GaAs. [The comparison of the model parabolic 
potential with the calculated potential in a realistic split 
gate wire is shown in Fig. 1(b)]. By varying Vq and uiq 
we can change the wire width and the electron density; 
in our calculations we set the Fermi energy Ep = 0. 

In order to calculate the magnetoconductance of the 
quantum wire, its subband structure and the wave- 
functions we use the Green's function techniqu o^^i^° 
where the electron interaction and the spin effects are 
included self-consistently within the framework of the 
Kohn-Sham density function theory in the local spin den- 
sity approximation^^. (The reliability of the the spin den- 
sity functional theory for electronic structure and mag- 
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netotransport calculations in quantum wires, dots and 
related structures is discussed in details in Refs. [23ll29h . 

We start from the Hamiltonian H = 
X^cr (-f^o + ^^'^(y)) , where Hq is the kinetic energy 
in the Landau gauge, 




and the total confining potential ^'^(y), 

V%y) = Feo„/(2/) + VH{y) + + (3) 

includes the bare electrostatic confinement Vconf{y) 
(given by Eq. ([1])), the Hartree potential VH'(y), the 
exchange-correlation potential V^^{y), and the Zeeman 
term, ~ gfibBa, where cr = ±i describes the spin- 
up and spin-down states, "f ,],; fif, ~ eh/2me is the Bohr 
magneton, and the bulk g factor of GaAs is g = —0.44. 
The Hartree potential Vff (y) due to the electron density 
'^(y) ~ 12(7 "■'^(y) (including the mirror charges) reads^^ 



Vniy) 



dy'n{y') In 



{y-v'Y 



{y - y'Y + 462 



, (4) 



with b being the distance from the electron gas to the 
surface (we choose h = 200 nm). The exchange and cor- 
relation potential Vxc{y) in the local spin density approx- 
imation is given by the functional derivative 



ySM^^{ne.c{n,ay))} 



(5) 



where C,{y) 



is the local spin-polarization. In 



the present paper we use the parameterization of the ex- 
change and correlation energy e^c given by Tanatar and 
Ceperly (TC)2£. Note that we also performed calcula- 
tions on the basis of the parametrization recently pro- 
vided by Attaccalite et al^ and found only marginal 
difference with the results based on the TC functional. 

The spin-resolved electron density in the wire can be 
expressed via the Green's function G'^{y,y,E) 



'{y) = --im 

TT 



dEG^iy,y,E)fFD{E-EF] 



(6) 

where fpoiE — Ep) is the Fermi-Dirac distribution func- 
tion. The Green's function, the Bloch states, the electron 
and current densities are calculated self-consistently us- 
ing the technique described in detail in Ref. [l^ . Knowl- 
edge of the wave vectors for different Bloch states a 
allows us to recover the subband structure, i.e. to cal- 
culate an overage position y'^ of the wave functions for 
different modes a for the given energy 

= — • (7) 

We calculate the spin-resolved conductance of the wire 
on the basis of the linear-response Landauer formula. 



dE 



dfjE-Ep) 
dE 



(8) 



where summation is performed over all propagating 
modes a for the spin a, with E^f^^ being the propaga- 
tion threshold for a-th mode. The current density for a 
mode a is calculated a.si^ 



dE- 



jl{y,E) ( df {E - Ef) 
dE 



(9) 



with and j^{y,E) being respectively the group veloc- 
ity and the quantum-mechanical particle current density 
for the state a at the energy E, and V being the applied 
voltage. All the calculations presented in this paper are 
performed for the temperature T = 100 mK. In order 
to speed up the calculation we use the modified Broy- 
den method^ that allows one to reduce the number of 
iterations need to achieve a self-consistent solution from 
- 2000 to only - 50 . 



III. RESULT AND DISCUSSIONS 

A. Hartree and spin DFT approximations 

We start our analysis of the magnetoconductance and 
the magnetosubband structure in quantum wires from 
the case of the Hartree approximation when the exchange 
and correlation interactions are not included in the ef- 
fective potential (i.e. when V^^{y) is set to zero in Eq. 
([3])). Note that the total potential V^ly) also includes 
the Zeeman term leading to the spin splitting even in the 
Hartree case. The effect of the Zeeman term is however 
negligibly small in the considered field intervals. We will 
thus refer to the Hartree case as for the case of spinless 
electrons. The results obtained in the Hartree approxi- 
mation will provide a basis for understanding of the effect 
of the exchange and correlation within the spin DFT ap- 
proximation. 

Figure [2)^a),(b) shows the magnetoconductance of a 
representative wire with the effective width w « 350 
nm and the electron density in the center of the wire 
n(0) = 3.2 • 10^^ m~2 calculated within the Hartree and 
the spin DFT approximations. The Hartree magnetocon- 
ductance shows the plateaus quantized in units of 2e'^/h 
separated by transition regions whose width grows as the 
magnetic field is increased. For large fields the width of 
the transition regions is comparable or can even exceed 
the width of the neighboring plateaus. For low fields, 
B < Bcrit, the width of the transition regions practically 
shrinks to zero; for the quantum wire at hand this critical 
field is i?crit ^ 0.6T, corresponding to the subband index 
~ 17, see Fig. [UJa). Note that in a standard one- 
electron picture of edge states the magnetoconductance 
of a clean wire (without impurities) is strictly quantized 
in units of 2e'^ /h (for spinless electrons), with vanishing 
width of the transition regions. Formation of the transi- 
tion region between the plateaus is shown to be related 
to development of the compressible strip in the middle of 
the wire^. 
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FIG. 2: (a) Conductance of the quantum wire calculated within the spin DFT and in the Hartree approximation (the later is 
shifted by —2e^/h for clarity). The vertical lines are drawn to emphasize that the width of the odd conductance steps in the spin 
DFT calculations is equal to the width of the transition intervals between the conductance steps in the Hartree calculations, (b) 
Comparison of between the Hartree magnetoconductance and the magnetoconductance calculated according to the Chklovskii 
et al. conjecture Gch = ^^v{Q). (c) Evolution of the magnetosubband structure in the interval 6e^ /h < G < 8e^ /h (second 
and third row). Fat solid lines indicate the total confining potential. The first row shows the electron density profiles (local 
filling factors) i^(j/) = n{y)/nB {ns = eB/h). Two lower rows show the current densities for the last two subbands {N = 7, 8). 
Parameters of the wire are the same as in Fig. [T] T = lOOmK. 
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FIG. 3: Solid line: The width of the transition regions be- 
tween the Hartree plateaus (which is equal to the width of 
the odd plateaus in the spin DFT calculations) as a function 
of 1/A*', with A'' being the subband number (indicated in the 
plot). Thin line: The subband splitting in the wire center 
AS = E (the definition of A^E is iUustrated in Fig. g}. Pa- 
rameters of the wire are the same as in Fig. [1] T = lOOmK. 



Let us now turn to the spin-resolved magnetoconduc- 
tance calculated by the spin DFT. The most striking fea- 
ture of the wire magnetoconductance is that the width of 
the odd conductance steps in the spin DFT calculations 
is equal to the width of the transition intervals between 
the conductance steps in the Hartree calculations, see Fig. 
dja). We will demonstrate below that the characteristic 
features in the spin-resolved conductance of the quantum 
wires calculated on the basis of the spin DFT (including 
the dependence of the width of the odd plateaus on the 
magnetic field and collapse of the odd plateaus at lower 
fields B < i?crit) can be understood from the analysis of 
the compressible strip structure for spinlcss electrons and 
the corresponding magnetoconductance and the magne- 
tosubband structure evolution in the Hartree approxima- 
tion. 

Before we proceed to the analysis of the evolution of 
the magnetosubband structure, it is instrumental to out- 
line how the exchange interaction induces the spin split- 
ting (see for details Refs. [l9ll20[ ). The spin splitting is 
most pronounced in the compressible strips. Indeed, the 
Hartree compressible strips are formed for partially occu- 
pied states in the vicinity of Ep when the Fermi-Dirac oc- 
cupation fp£i < 1 (i.e. in the window \E — Ef\ < 2'KkT). 
When the states are partially occupied, the system be- 
haves like a metal, where the electrons can easily read- 
just their density to screen the external potential. It is 
important to stress that the compressible strips, being 
partially occupied, allow for different population of the 
spin-up and spin-down states. In the DFT calculations, 
this population difference (triggered by the Zeeman split- 
ting) is strongly enhanced by the exchange interaction. 
This leads to the lifting of the subband degeneracy and 
to the spatial separation between the spin-up and spin- 
down states. 



B. The magnetosubband structure and the 
magnetoconductance 

In order to get insight into evolution of the odd con- 
ductance plateaus let us inspect the Hartree and the spin 
DFT magnetosubband structure. Let us, for example, 
concentrate at the field region where 6e^/h < G < Se'^/h, 
i.e. when the magnetoconductance clearly shows the spin 
splitting. The magnetosubband structure for several rep- 
resentative fields in this region is shown in Fig. [DJc). In 
all subsequent discussions we will focus on the two high- 
est subbands (in this case = 8 and N = 7), because 
the depopulation of these two subbands determines the 
features in the conductance steps (note that all the re- 
maining subbands are fully filled). At B ~ 1.62T the 
subbands N = 7,8 are fully occupied and thus the to- 
tal conductance G = 8e'^/h. The Hartree calculations 
for spinless electrons show the presence of a narrow com- 
pressible strip of the width w^^p, see Fig. [5] A. When 
the exchange interaction is included, the subbands split 
which leads to the spatial separation between the spin-up 
and spin-down states. (The spatial spin separation due 
to the suppression of the Hartree compressible strips was 
discussed in detail in Ref. [20[ ). 

When magnetic field is increased the subbands are 
pushed up in energy (see Fig. [2]B; i? = 1.7T). The com- 
pressible strip in the Hartree calculations becomes wider 
(because the confinement is smoother in the wire center) , 
and it moves closer to the center of the wire. The ex- 
change interaction quenches the compressible strip caus- 
ing the splitting of the spin-up and spin-down subbands. 
However, despite of the lifting of the spin degeneracy the 
subband bottoms arc still below Ep at the wire center. 
Because of this, the two spin split subbands iV = 7, 8 re- 
main fully (and equally) populated and the conductance 
remains on the plateau G ~ 8e^//i. 

When the magnetic field is increased to _B = 1.78T 
(Fig. mC), the Hartree compressible strip reaches the 
middle of the wire. This means that the subbands be- 
come partially occupied because their bottoms are now 
within the window \E — Ep\ < 2'KkT (where fpD < 1). 
As a result, the conductance of the spinless Hartree elec- 
trons starts to decrease and the transition region between 
the plateaus starts to form. With further increase of 
the magnetic field the Hartree compressible strip in the 
middle of the wire shrinks, and at S 2.05T the sub- 
bands depopulate completely, as they are pushed above 
the window \E — Ep] ~ 27rfcT. We conclude the dis- 
cussion of the evolution of the Hartree subbands by re- 
emphasizing the fact that the transition between the con- 
ductance steps starts when the compressible strip reaches 
the center of the wire and it ends when the compress- 
ible strip disappears and two highest (spin-degenerate) 
magnetosubband are pushed above Ep. Note that even 
though this discussion was focused on the transition be- 
tween G = 8e'^/h and 6e^//i-plateaus in the Hartree con- 
ductance, the same scenario of the Hartree subband de- 
population holds for all other subbands. 
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Let us now examine how the exchange interaction af- 
fects the transition region between the Hartree plateaus 
G = 8e^//i and Qe^ /h. Similarly to the cases of lower 
fields discussed above (Fig. S] A, B), the exchange in- 
teraction causes the subband repulsion and the spatial 
spin separation of the wavefunctions (the latter being 
equal to the width of the Hartree compressible strips, 
see the lower panel of Fig. Wc)). For B = 1.78T the 
Hartree compressible strip covers the central part of the 
wire. As a result, the bottom of the higher energy (spin- 
down N = 8) subband is situated within the window 
\E — Ef\ < 27rfcT (and thus this subband is only par- 
tially populated) , whereas N = 7 (spin-down) subband is 
pushed below Ep and thus remains fully populated (Fig. 
4 C, spin DFT calculations). Thus, at B = 1.78T the 
transition to the odd plateau G = 7e^//i starts to form. 
The exchange interaction keeps 8th and 7th subbands 
separated such that with further increase of the magnetic 
field 8th subband becomes quickly depopulated while 7th 
subband remains fully occupied with its bottom being be- 
low Ep, see Fig. 4 D, E (B = 1.84T and B = 2.05T). 
The above field interval (i.e. 1.84T < B < 2.05T) corre- 
sponds to the odd step in the magnetoconductance. Fi- 
nally, at B = 2.05T (i.e. at the same field when the cor- 
responding Hartree subbands depopulate), the bottom 
of 7th subband is pushed above Ep (to be more precise, 
above Ep + inkT) and 7th plateau in the conductance 
disappears. 

To summarize the discussion presented in this section, 
the formation of the odd magnetoconductance plateaus 
due to the exchange interaction can be traced to the for- 
mation of the compressible strips in the center of the 
wire in the case of the spinless electrons. The exchange 
interaction lifts the spin degeneracy such that the bot- 
tom of the highest (even) subband remains pinned to Ep^ 
whereas the bottom of the highest odd subband remains 
below Ep- As a result, the odd plateaus (whose width is 
equal to the width of the transition regions between the 
Hartree plateaus) develop in the magnetoconductance. 

Note that the analytical solution to the electrostatic 
problem of the electron density distribution in a quan- 
tum wire for spinless electrons has been obtained by 
Chklovskii, Matveev and Shklovski^ for the high mag- 
netic field regime (when only a few lower Landau levels 
are occupied, iV < 2 ~ 4). [A good agreement with with 
the analytical results of Chklovskii et al. has been re- 
ported by Oh and Gerhardts within the self-consistent 
Thomas- Fermi calculations^]. Chklovskii et al^ have 
also discussed the magnetoconductance of the quantum 
wire. They found that in a realistic quantum wire the 
conductance plateaus are practically absent, i.e. the 
conductance is not quantized (see Fig. 5 in Ref. 0). 
This conclusion is in obvious disagreement with the ex- 
perimental results showing pronounced plateaus in the 
two-terminal magnetoconductance at integer values of 
e^//i.— 1^1^. They attributed this discrepancy to the pres- 
ence of disorder in the channel. We however, have demon- 
strated above that even in an ideal clean channel (without 



disorder) the conductance shows the pronounced quan- 
tization with wide plateaus and sharp rises. The reason 
of the discrepancy we instead attribute to the conjecture 
used by Chklovskief al. that the ballistic conductance 
is given by the filling factor in the middle of the wire, 

Gch ~ ^IT^^^)- ^^'^ quantum mechanical calculations 
show that this conjecture is not justified. Indeed, the 
comparison of the conductance calculated according to 
the Chklovskii et al. conjecture with the magnetocon- 
ductance calculated according to Eq. ^ (see Fig. El^b)) 
shows that Gch reproduces the overall decrease of the 
conductance rather well, but does not at all recover the 
steps in the conductance related to the subband depop- 
ulation (see also Ref. [3^ for a related discussion). Our 
results thus indicate that while electrostatic and Thomas- 
Fermi-type approaches can be very successful in the de- 
scription of the electron density and the structure of the 
compressible and incompressible strips for the spinless 
electrons, an accurate description of the magnetoconduc- 
tance requires detailed quantum mechanical information 
for the wave functions and the currents densities. 



C. Collapse of the odd magnetoconductance 
plateaus at lower fields 

When magnetic field is increased, the width of the 
transition regions between the Hartree plateaus (which 
is equal to the width of the odd plateaus in the spin 
resolved magnetoconductance), AS, is also gradually in- 
creases, see Figs l2l3l We attribute this increase of /S.B to 
the effect of the enhanced electron screening due to the 
evolution of the compressible strip in the middle of the 
wire. Indeed, the transition regions between the Hartree 
plateaus are related to the depopulation of the high- 
est (spin-degenerate) subbands forming the compressible 
strip in the center. In high magnetic field each subband 
(representing a Landau level) accommodates the same 
number of electrons, such that the density of the elec- 
trons in the highest subband is proportional to ~ 1/N. 
Thus, one can expect that the width of the compressible 
strip in the middle of the wire Wcomp and hence the width 
Ai? grow as B increases (note that B ^ ^/N). Figures[3] 
- [5] illustrating the magnetic field dependence of AS and 
Wcomp confirms this expectation. Note that AB shows 
a nonlinear dependence on That is, for low fields, 

B < Bait, the width AB rapidly decreases when B de- 
creases, such that the odd plateaus are no longer seen in 
the magnetoconductance. 

Let us now concentrate on this feature of magnetocon- 
ductance in more detail. Figurc[2]shows the spin-resolved 
magnetoconductance G^ and G^ . It is worth to stress 
that the spin degeneracy remains lifted even for fields 
smaller than Bcrit ■ The total conductance G = G^ + G^ 
does not however exhibits the odd plateaus for B < Bcrit 
because the strength of the exchange splitting becomes 
comparable to the thermal broadening of the plateaus. 
This is illustrated in Fig. [3] that shows the dependence 
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FIG. 4: The magnetosubband structure of a quantum wire 
within the Hartree and the spin-DFT approximation for dif- 
ferent number of occupied subbands A'^. Fat sohd lines indi- 
cate the total confining potential. [For all cases the magnetic 
field is chosen such that the Hartree compressible strips in 
the middle of the wire have a maximal width] . AE shows the 
subband splitting in the center of the wire. Parameters of the 
wire are the same as in Fig. [1] T = lOOmK. 



of the subband splitting AE in the center of the wire 
on the magnetic field and its comparison to the energy 
window AirkT (where the derivative of the Fermi-Dirac 
distribution function is distinct from zero). [The defini- 
tion of the subband splitting AE is outlined in Fig. 0]. It 
is also worth pointing out that in accordance to the pre- 
vious discussion, AE and AB exhibit similar behavior as 
a function of magnetic field, see Fig. [S] 

In order to understand the nonlinear behavior of AB 
leading to quenching of the odd magnetoconductance 
plateaus at low field we examine the wave functions and 
the current density distributions. Figure [5] (a) shows the 
current density for the Hartree subbands = 8 — 18 
along with the maximal width of the Hartree compress- 
ible strip Wcomp in the middle of the wire. The ex- 
tend of the wave function for the highest iVth subband, 
{iPn) VnIb, [Ib is the magnetic length) gradually 
increases when the magnetic field is lowered^. Note 
that (V'jv) is larger than the width of the compressible 
strip Wcomp already for = 8. When the extend of 
the wave function exceeds the width of the compressible 
strip, the ability of the system to screen the external po- 
tential is greatly reduced because the wave function can 
be shifted within the distance not exceeding the width 
of the compressible strip Wcomp- Thus, the smaller the 
ratio Wcomp/ {''Pn) is, the weaker is the effect of the redis- 
tribution of the electron density required to screen the 
external potential. This reduced screening efficiency for 
lower fields (when WcompI (V'Af) 1) translates into the 
suppressed exchange splitting and thus to disappearance 
of the odd magnetoconductance plateaus. 
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FIG. 5: (a) The current densities for the highest occupied 
Hartree subband compared to the maximal width of the com- 
pressible strip in the middle of the wire for different magnetic 
fields. A'' indicates the subband number (see Fig. |4] for the 
corresponding magnetosubband structure). (b)The current 
densities for the highest occupied Hartree subband (A'' = 8) 
compared to the maximal width of the compressible strip in 
the wires with the effective widths of u; = 350 nm and w — 700 
nm. Parameters of the wire are the same as in Fig. [1] param- 
eters of the second wire are indicated in Fig. [S] T = lOOmK. 



Note that the extend of the wave function (iPn) for a 
given subband number N (or for a given magnetic field) 
is not particularly sensitive to the wire width w (at least 
in the regime when the cyclotron radius rc < w). At the 
same time, the maximum width of the compressible strip 
increases with increase of the wire width. This is illus- 
trated in Fig. Elb) that shows the current density dis- 
tribution and Wcomp for two quantum wires of the with 
w = 350 nm and 700 nm for the case of A^ = 6 occu- 
pied subbands. [Note that in the bulk limit (i.e. for 
the edge of the 2DEG) the compressible strip covers the 
semi-infinite space, such that regardless of the subband 
number, Wcomp/ i'ipN) 3> !]• Therefore, for the given N 
(or magnetic field B) the ratio Wcomp/ (^w) is larger in a 
wider wire and therefore the screening efficiency is higher. 
One can therefore expect that in a wider wire the mag- 
netosubband spin splitting due to exchange interaction 
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FIG. 6: Comparison between the calculated magnetoconduc- 
tance of the wire with the effective width of it; = 70Qnm and 
the experimental magnetoconductance^. T = lOOmK. 



leading to the appearance of the odd magnetoconduc- 
tance plateaus would manifest itself for larger subband 
numbers (lower fields). Our calculations show that this 
is indeed the case. For example, in a wire with w = 350 
nm the odd plateaus become discernible for the subband 
index N = 17 (i?crit ^ 0.6T), whereas for the wire with 
w = 700 nm the last odd plateau is seen for N ^ 19 — 21 
(Bent 0.4T), c.f. Figs. [Handin] 



IV. COMPARISON TO THE EXPERIMENT 

We have performed magnetotransport calculations for 
several quantum wires with effective widths in the range 
200-700 nm and the electron densities 1.5 — 3.5- lO^^m"^. 
All the wires exhibit the same behavior described in de- 
tail in previous section. 

A detailed comparison of our calculations with the ex- 
perimental magnetoconductance^ for a some representa- 
tive quantum wire is shown in Fig. [SI The width of 
the wire is estimated to be ~ 680 nm, and the sheet 
electron density in the bulk Us = 2.15 • 10^^m~^fi. The 
theoretical magnetoconductance shown in Fig. [6] is cal- 
culated for the electrostatic confining potential with the 
parameters Vq = —0.4 eV, Huiq — 1.91 meV (giving the 
effective wire width w = 700 nm and the electron density 
in the wire center n = 2.2 • lO^^m^^). In the both ex- 
perimental and simulated wires the electrons is situated 
at the distance ~ 200 nm below the surface. The com- 
parison to the subband depopulation in the experimental 
structure demonstrates that such a choice of the param- 
eters provides a satisfactory approximation for the ac- 



tual confining potential. We stress here that a magnetic 
field dependence of the subband depopulation can be de- 
scribed by the one-particle Schrodinger equation (for a 
given confining potential)^, whereas our main focus here 
is the electron interaction effects leading to formation 
of the odd steps in the magnetoconductance due to the 
exchange interaction. The comparison of the calculated 
and the experimental curves demonstrates a good quan- 
titative agreement between the widths of the odd (as well 
as even) plateaus in the calculated and the experimen- 
tal magnetoconductance. The calculations also provide a 
reasonably close estimation of the subband index corre- 
sponding to the last resolved odd plateau in the magne- 
toconductance, ~ 19 — 21, whereas the corresponding 
experimental value is A^ = 15. 

It would be unreasonable to expect an exact agreement 
between the theory and the experiment. There are sev- 
eral factors that have not been taken into account in the 
theoretical modelling. We list below some of them. 

(a) The experiment^ is performed in the QPC geome- 
try whereas our calculations are done for an infinite quan- 
tum wire. In the edge state transport regime considered 
here with A^ 3> 1 this is not expected to be a source 
of significant discrepancy between the theoretical predic- 
tions and the experiment. Nevertheless, the effect of the 
QPC geometry on the magnetoconductance remains to 
be seen. 

(b) The calculations are performed for ideal clean 
wires without impurities. In a one electron descrip- 
tion the transition regions between the magnetoconduc- 
tance plateaus are the step functions with a zero width. 
A random impurity potential is known to lead to the 
smooth transition regions of a finite width even in the 
one-electron picture^. A comparison of the magneto- 
conductance traces in Fig. [H] clearly shows that the the 
transition regions seen in the experiment are significantly 
wider than the theoretical ones. [We remind that a finite 
width of the transition regions in the theoretical magne- 
toconductance is due to the formation of the compressible 
strips in the middle of the wire]. We attribute the dif- 
ference in the widths of the calculated and experimental 
transition regions to the effect of the disorder that has 
not been included in the model. Note the experimen- 
tal magnetoconductance traces for narrower QPCs show 
the conductance fluctuations in the transition regions^, 
which is a clear manifestation of the disorder potential 
due to impurities^. Note that the presence of disorder 
can also lead to the destruction of the exchange enhance- 
ment of the g-factor and thus to the collapse of the spin 
splitting (i.e. to the suppression of the odd plateaus for 
A^ > Nc)^. This effect does not seem to be relevant 
to the experiment because the spin splitting, being sup- 
pressed in the narrow structures (QPCs) for larger A^, is 
still clearly seen in the bulk Hall measurements. 

(c) In our calculations we assumed that an elec- 
tron motion is confined to a two-dimensional plane, 
which is a good approximation for heterostructures where 
the electrons are localized on the interface between 
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GaAs/AlGaAs. In the experimental structures^ the elec- 
trons are confined in a quantum well that is populated by 
donors situated on both sides from the well. An accurate 
description of this geometry might require accounting in 
the Schrodinger equation the electron motion in the di- 
rection perpendicular to the interface. 

Finally, we notice that while we compared our calcu- 
lations with the experimental conductance of one repre- 
sentative wire, our spin-DFT calculations qualitatively 
reproduce all the features observed in other samples. 
This includes the dependence of of the width of the odd 
plateaus AS on the magnetic field shown in Fig. [31 It 
is also worth to stress that the theory confirms (and ex- 
plains) the experimental finding that in wider wires the 
collapse of the odd plateaus occurs at lower fields. We 
however are not in position to fit all the experimental 
data. This is because that this task would require a de- 
tailed knowledge of the bare confining electrostatic po- 
tential Vconf due to the gates and the donor layers, Eq. 
([I]). Indeed, the bare electrostatic confinement deter- 
mines the total self-consistent confining potential (y) , 
which, in turn, determines the depopulation of the mag- 
netosubbands (i.e. the dependence of the subband num- 
ber N on B)^'^. We are not in a position to perform 
a systematic search for the parameters of Vconf giving 
rise to the B-dependence of the subband depopulation 
consistent with each experimental magnetoconductance 
trace. This is simply because of a computation burden 
related to this task: each point on the magnetoconduc- 
tance plot requires up to one hour of a processor time. 
[Note that a calculation of the electrostatic confinement 
and the self-consistent potential starting from the layout 
of the actual hctcrostructurc of Ref. represents a sepa- 
rate task which is outside a scope of the present study] . 



V. CONCLUSION 

In this paper we provide a systematic quantitative de- 
scription of the magnetoconductance of the split-gate 
quantum wires focusing on the formation and evolution 
of the odd conductance plateaus. In order to calculate 
the electron density, magnetosubband structure and the 
magnetoconductance we utilize the self-consistent Greens 
function technique combined with the spin density func- 
tional theorjii^. 

We start our analysis with the case of spinless elec- 
trons in the Hartree approximation (disregarding the 
exchange and correlation interactions). The calculated 
Hartree magnetoconductance shows the plateaus quan- 
tized in units of 2e^//i separated by transition regions 
whose width grows as the magnetic field is increased. The 
transition regions arc attributed to the formation of the 
compressible strips in the middle of the wire occupied 
by electrons belonging to the highest (spin-degenerate) 
subband. In agreement with experiments, the width of 
the transition regions for large fields is comparable to the 
width of the neighboring plateaus. This is in contrast to 



both the one-electron description where the conductance 
shows the step-like behavior with the rises between the 
plateaus of zero width, as well as to the electrostatic the- 
ory of Chklovskii et al^ where the magnetoconductance 
exhibits narrow plateaus of negligible width separated by 
much broader transition regions where the conductance 
is not quantized. 

Accounting for the exchange and correlation interac- 
tions within the spin DFT leads to the lifting of the spin 
degeneracy and formation of the spin-resolved plateaus 
at odd values of /h. The most striking feature of the 
magnetoconductance is that the width of the odd con- 
ductance steps in the spin DFT calculation is equal to 
the width of the transition intervals between the conduc- 
tance steps in the Hartree calculations. This is because 
the transition intervals in the Hartree magnetoconduc- 
tance correspond to the formation of the compressible 
strip in the middle of the wire. At the same time, in 
the compressible strip in the center of the wire the states 
arc only partially occupied. As a result, the exchange 
interaction enhances the difference in the spin-up and 
spin-down population, which leads to the lifting of the 
subband spin degeneracy and formation of the odd con- 
ductance plateaus. 

In agrcmcnt with the experimental results^, we find 
that the width of the odd magnetoconductance plateaus 
gradually decreases with decrease of the magnetic field. 
For lower fields B < Bent, the odd plateaus rapidly disap- 
pear such that the magnetoconductance shows the quan- 
tization in units of 2e^/h. The wider the wire, the lower 
the critical field Scrit corresponding to the disappearance 
of the last resolved odd plateau. We attribute this effect 
to the reduced screening efficiency in the confined (wire) 
geometry when the width of the compressible strip in 
the center becomes much smaller than the extend of the 
wave function. This, in turn, leads to the suppressed 
exchange splitting and collapse of the odd magnetocon- 
ductance steps. 

A detailed comparison to the experimental datai (see 
Fig. [6]) demonstrates that the spin-DFT calculations 
reproduce not only qualitatively, but rather quantita- 
tively all the features observed in experiment. This in- 
cludes the dependence of the width of the odd and even 
plateaus on the magnetic field as well as the estimation 
of the subband index corresponding to the last resolved 
odd plateau in the magnetoconductance. The experi- 
ment however shows wider rises between the transitions 
plateaus in comparison to the calculated ones. We at- 
tribute this difference to the effect of smooth potential 
due to remote donors that has not been accounted for in 
our calculations (performed for clean disorder-free wires). 
Despite of this discrepancy, the overall good agreement 
between the theory and experiment makes it possible to 
conclude that the spin DFT approach represents the pow- 
erful tool to study large realistic quantum Hall systems 
containing hundreds or thousands of electrons, providing 
detailed and reliable microscopic information on wave- 
functions, electron densities and currents as well as the 
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conductance. 
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